Extended quantum critical phase in a magnetized spin-| antiferromagnetic chain 
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Measurements are reported of the magnetic field dependence of excitations in the quantum crit- 
ical state of the spin S = 1/2 linear chain Heisenberg antiferromagnet copper pyrazine dinitrate 
(CuPzN). The complete spectrum was measured at fesT/J < 0.025 for H = and H = 8.7 Tesla 
where the system is ~ 30% magnetized. At H = 0, the results are in quantitative agreement 
with exact calculations of the dynamic spin correlation function for a two-spinon continuum. At 
high magnetic field, there are multiple overlapping continua with incommensurate soft modes. The 
boundaries of these continua confirm long-standing predictions, and the intensities are consistent 
with exact diagonalization and Bethe Ansatz calculations. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.50.Ee 



One of the most important ideas to emerge from stud- 
ies of condensed matter systems in recent years is the 
concept of quantum criticality A quantum critical 
point marks a zero temperature phase transition between 
different ground states of a many-body system as a re- 
sult of changes in parameters of the underlying Hamilto- 
nian. Precisely at the quantum critical point the system 
is without characteristic length scales or energy scales, 
with power-law spatial correlations and gapless excita- 
tions. Finite temperature properties close to quantum 
criticality are anomalous and reflect universal properties 
of the underlying quantum field theory. 

While fine tuning of a parameter in a system's Hamil- 
tonian is generally required to achieve quantum critical- 
ity, it is inherent to the spin S=l/2 linear chain Heisen- 
berg antiferromagnet (LCHAFM). The spin dynamics of 
the LCHAFM have been studied in a number of mate- 
rials |2|, y, U, [5J, |fj , and the elementary excitations are 
S=l/2 spinons that form a gapless, two-particle contin- 
uum 0>13- In a magnetic field, H the Hamiltonian of 
the LCHAFM is H= ^ t [./S,S l+J - g f i B HSf}. While 
changes in Ji such as the introduction of dimerization 
or interchain coupling drive the LCHAFM away from 
criticality |9(, the system should remain quantum crit- 
ical at T = in fields below the full y s pin polarized 
state, which occurs at H c = 2 J/ gnu '.10, 11]. Promi- 
nent among the features predicted to exist along this 
quantum critical line are a set of field-dependent two- 
particle continuua 0, 0, E3, 0] with incommensurate 
soft modes that move across the ID Brillouin zone with 
increasing field as = 2irm and (7^2 = 7r — 27tto, where 
< to < 1/2 is the reduced magnetization per spin. 

Experiments on copper benzoate [l5| have verified 
the predicted field-dependence of the incommensurate 
wavevector 5^2 • However, due to a combination of a stag- 
gered g-tensor and Dzayaloshinskii-Moria interactions, 



the field drives that system away from the critical line 
to a state with confined spinons and a gap in the exci- 
tation spectrum |l6| . In contrast, the finite-field criti- 
cal state of the S=l/2 LCHAFM is accessible in copper 
pyrazine dinitrate, Cu(C 4 H 4 N2)(N03)2 (CuPzN). This 
well characterized organo-metallic magnet has J = 0.9 
meV, and negligible interchain coupling (J'/ J < 10~ 4 ) 
|6ul7|. The spin chains in CuPzN are uniform, with one 
Cu 2+ ion per unit cell along the chain, and specific heat 
measurements have shown that CuPzN remains gapless 
for H < Q.6Hc — 9 T Q. In this paper we report inelas- 
tic neutron scattering measurements of the full spectrum 
of CuPzN, both at H = and at H = 8.7 T, where 
771 = 0.15. In zero field, the spectrum is consistent with 
the exact two-spinon contribution to the spin fluctuation 
spectrum [||. The finite field data show the long-sought 
field-dependent continua 7], and are in detailed agree- 
ment with theoretical and numerical work 
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The measurements were performed using the SPINS 
cold neutron triple axis spectrometer at the NIST Cen- 
ter for Neutron Research. The sample studied contained 
22 single crystals of deuterated CuPzN ^3 with a total 
mass of 3.06 grams, and a measured mosaic of two de- 
grees. We used a dispersive analyzer configuration that 
detects scattered neutrons with an angular ran ge o f 10° 
and an energy range 2.7 meV < Ef < 3.55 meV |l9j. We 
used a Be filter before the sample for < 5.15 meV, and 
a BeO filter after the sample. The measured full width 
at half maximum (FWHM) elastic energy resolution was 
Shu> — 0.14 meV, and the FWHM wave-vector resolu- 
tion along the [100] chain axis was 5Q|| = 0.03 A" 1 at 
Q = (hkl) = (fjjO) at hoj = [1J. Data were obtained 
along the line between Q = (l|0) and Q = (2^0), with 
H\\c. 

Figure OJa) shows the normalized magnetic scattering 
intensity X m (q,u>) for CuPzN, measured in zero field at 
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FIG. 1: (a) Magnetic inelastic neutron scattering intensity 
T m for CuPzN at T = 0.25 K and H = versus wave- vector 
and energy transfer, (b) Calculated two-spinon contribution 
to X m . Solid lines are the predicted lower and upper bounds 
of the spinon continuum. A representative FWHM resolution 
ellipsoid is shown in (a). 



T = 0.25 K. Wave-vector transfer along the chain is rep- 
resented as q = 2w(h — 1), This data set was obtained by 
combining data taken at Ei = 3.35 meV, 3.55 meV, and 
for 3.75 meV < E z < 5.75 meV with 40 fiV steps. After 
subtracting the background scattering measured with the 
analyzer in a non-reflecting geometry, the incoherent elas- 
tic scattering profile of the sample was determined from 
data that were at least 0.2 meV outside of the known 
bounds || of the H = spinon continuum. This profile 
was scaled to to the measured elastic incoherent scat- 
tering intensity at each q, and subtracted from the raw 
data. The data were converted to the normalized scat- 
tering intensity by comparison with the measured inco- 
herent scattering intensity of a vanadium standard. The 
data at 2n — q were then averaged with that at q, binned 
in bins of size dhuj = 25 /xeV by 5q = 0.026-7T, and av- 
eraged over a rectangle of the same size as the FWHM 
energy and wave-vector resolutions to produce the color 
contour plot shown in Fig. }Ua). 

These data provide a complete picture of the zero-field 
spinon continuum (SC) in CuPzN. The solid lines are 
the SC's lower and upper bounds, calculated for J = 0.9 
meV |7j. Figures Eta)-(c) and[2[a)-(c) show cuts through 
the data, and highlight several important features. These 
include the monotonically decreasing intensity at q = tt 
[Fig. Il^a)], the peaks near q = ir associated with the di- 
vergence in S(q, to) at the SC lower bound [j| [Fig|s2a)- 
(c)], and the smaller feature at low q where the SC nar- 
rows as q — > 0. Figures |2Ib)-(c) show the asymmetric 



H = H - 8.7 T 




12 12 

fico (meV) (meV) 



FIG. 2: Magnetic scattering intensity vs Tkj at constant q for 
CuPzN at T = 0.25 K and H = (a)-(c), and at H = 8.7 
T (d)-(f). Each point includes scattering within a width of 
Aq = 0.05-7T. The solid (dashed) lines in (a)-(c) are the exact 
(approximate) model for S(q,u) described in the text. 

lineshapes produced by the SC for q ^ tt. Figure |3fc) 
also shows data taken at T = 35 K. The decrease in in- 
tensity at q = ir confirms that the low-T scattering is 
magnetic. 

Figure^b) shows T m {q,oj) calculated 0] from the re- 
cently derived exact two-spinon contribution to S{q, lo) 
HHillll with J = 0.9 meV, with the addition of a con- 
stant intensity fixed at the average residual background 
measured 0.5 meV away from the continuum. X m (q,u>) 
is also shown as solid lines on the cuts in Figs. Efa)-(c) 
and[2a)-(c). There are no adjustable parameters in the 
model, and the agreement with the data is excellent, al- 
though we note that there is an overall 10% uncertainty 
in the vanadium normalization. 

Figures |21 and |21 also show a comparison (dashed lines) 
to the lineshapes produced by a global fit to the approx- 
imate form for S(q, lo) that has been previously used 
to model measurements of the H — SC in CuPzN Q 
and other materials 0, This approximation is now 
known to overestimate S(q, to) near the upper SC bound- 
ary 0,0]. Indeed, some indications of this effect are seen 
in Figs. Gib) and|2Ib), and although our data are not 
optimal for observing these differences, the exact model 
does give an overall better description of the data. 

The spectrum of CuPzN in a magnetic field is consider- 
ably more complex than at zero field, as seen in Fig. Ufa), 
which shows X m measured at H — 8.7 T. This data set 
combines data taken in the range 3.35 meV < Ei < 4.75 
meV with 20 fiV increments, and at Ei — 5.15 meV. An 
identical background subtraction, averaging, and binning 
procedure as used at H = was applied to these data, 
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FIG. 3: Magnetic scattering intensity vs q at constant Huj for 
CuPzN at T = 0.25 K at H = (a)-(c), and at H = 8.7 T (d)- 
(f). The data for q < n have been averaged with the data at 
2-K—q. The data displayed for q > n have not been averaged in 
this manner. Each point includes scattering within a width 
of Ahu = 0.1 meV. The solid (dashed) lines in (a)-(c) are 
the exact (approximate) model for S(q,u>) described in the 
text. The scattering measured at T = 35 K is shown as open 
symbols in panel (c). 



using the H = incoherent elastic profile. 

Several continua and prominent features have emerged 
that were not present at H — 0. These are highlighted 
in cuts shown in Figs. Efd)-(f) and |2Jd)-(f). First, a 
continuum of scattering is still observed at q = tt down to 
the lowest energy probed, which demonstrates that the 
system indeed remains gapless and critical at H = 8.7 
T. The spectrum has generally shifted to lower energy, 
and the strong ridge of scattering with linear dispersion 
near q = tt has a smaller slope than at H = 0, showing 
explicitly that the velocity of the elementary excitations 
has decreased [fjj. 

There is a strong peak in X m (g, u>) centered at q = tt at 
the field energy Tiw = gfi B H « 1.1 meV (g c = 2.07 
A weaker peak at this same energy can also be seen at q = 
0, which corresponds to uniform spin precession. Moving 
away from q = tt, there is another ridge of scattering 
intensity, which decreases in energy towards the expected 
field induced incommensurate wavevector q^2 = Q.Ttt. As 
shown in Fig. EJe), the incommensurate mode is seen in 
increased scattering at the lowest energies probed at 
compared to zero field [Fig.^b)] . Finally, there is a mode 
that begins at the field energy at q = 0, and decreases in 
energy with increasing q, but which loses intensity close 
to q~i,i = 0.37T, consistent with numerical work [12j. 

Using the Bethe Ansatz, Muller et al. identified six 
classes of excitations out of the partially magnetized 




FIG. 4: (a) Inelastic neutron scattering intensity T m (q,u>) for 
CuPzN at T = 0.25 K and H — 8.7 T. (b) Calculations of 
the different components of S{q,uj) for TV — 26 spins and 
m = 2/13 (H = 8.67 T for CuPzN). The area of each circle is 
proportional to S(q,u). (c) T m {q,ui) calculated for ensemble 
of chains with N = 24, 26, and 28. Lines show bounds of 
excitation continua for the finite- length S = 1/2 LCHAFM. 
Solid lines: continua predicted to predominate as N — » oo. 



ground state of an TV-spin chain that can contribute to 
w ) H- Each class produces a continuum, and Muller 
et al. determined their approximate boundaries. How- 



ever, three of these are predicted to dominate S(q,u)) as 
N — > oo, with (adopting the notation of Ref. |7J,) class 

(ii) contributing to S zz {q, u>), class (iii) to S u>), and 

class (vi) to 5_| (q, uS). The boundaries of these three 

continua, £2, £3, and £§, are shown as solid lines in Fig.0] 
They closely track the principal features of the data. 

We have calculated S(q, u>) for chains of length TV = 24, 
26, and 28, using the Lanczos technique G3. FigurelUb) 
depicts the results for N = 26 and m = 2/13, corre- 
sponding to H = 8.67 T for CuPzN 0. Each eigen- 
state of the chain is marked by a circle with an area 
proportional to the corresponding contribution to S(q, lo) 
UJ. The circles are color-coded to indicate which com- 
ponent of S(q,uS) the state contributes to. Note that 
in our experimental geometry we measure T m (q,u>) oc 

S Z z + ji^A + <5 h)- The corresponding continuum 

boundaries are also shown color-coded. Some spectral 
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weight for this finite chain is seen outside of continua 
(ii), (hi) and (vi), and so for reference the boundaries of 
the other three hnite--/V continua, £\, £^4, and £5, are in- 
cluded as dashed lines. The finite-chain results reproduce 
all trends and features of the measured intensity, and to- 
gether with the continuum bounds, suggest polarization 
assignments. 

A more direct comparison is achieved by combining the 
calculations in Fig.^b) with results for N — 24, m = 1/6 
and N = 28, m = 1/7, which correspond to H = 9.18 T 
and 8.21 T for CuPzN, respectively, and greatly increase 
the number of wavevectors sampled. The results at each 
iV were smoothed by convolving them with the response 
function of a finite-size system, so that 



S a/3 (q, huj) = 2J ^ S af 3 :N (qN, hu) 



N q N 



sm((q-q N )N/2) 



sm((q-q N )/2) 



, (1) 

The calculated results were then converted to I m (q,u>) 
0, and binned and averaged as described above. Af- 
ter multiplication by an overall scale factor, we obtain 
the contour plot shown in Fig. 0fc) which represents the 
measured intensity expected from an equal weighted en- 
semble of 24, 26, and 28 membered spin chains. While 
the relatively short chains yield stronger discrete energy 
bands than in the measurements, the simulation clearly 
captures the main features of the data. 

Karbach and Muller have recently identified a new 
class of quasiparticles for the partially magnetized S=l/2 
chain. These "psinons" play a similar role in the spec- 
trum as do the spinons at H = 0, yielding a continuum 
for each component of S(q,u>) [lijllj], with boundaries 
that are consistent with the approximate analytic ex- 
pressions that were used above. Karbach et al. have 
computed the psinon lineshapes for m — 0.25, where 
%i = &,2 = tt/2. This occurs at H = 11.9 T for CuPzN, 
but a qualitative comparison can still be made, as the 
properties of the lineshapes should vary smoothly. Of 
particular interest is the peak at the upper boundary of 
the psinon continuum for S zz (q, u>), which is due to a sin- 
gularity in the psinon density of states. This may explain 
the line of scattering intensity in Fig. Ua) that tracks the 
top of the £2 continuum. This is notably different from 
the zero-field case, where the corresponding singularity 
in the spinon density of states is compensated by a van- 
ishing matrix element, and S(q,u>) vanishes at the upper 
boundary of the spinon continuum, as seen in Fig. ^ 

Finally, we note that Fig. Ufa) shows some evidence 
of weak scattering intensity for Tilo > 2 meV. This could 
possibly be due to the presence of short chains resulting 
from impurities, or to higher-order processes not included 
in the spinon/psinon picture. However, we note that our 
error bars are much larger here than at lower energy due 
to shorter counting times, (see Fig.EJ, and so a definitive 
statement on the existence of excitations in this energy 
range cannot be made at this time. 



In summary, our experiments in CuPzN provide the 
first experimental example of an extended critical state 
in a quantum magnet. This detailed mapping of the spin 
excitation spectrum in the spin-1/2 linear chain antifer- 
romagnet verifies long-standing predictions based on the 
Bethe Ansatz of a field driven and critical incommensu- 
rate state. There is excellent agreement with finite chain 
calculations, and good qualitative agreement with the 
lineshapes predicted at higher fields based on novel "psi- 
non" quasiparticles. Recent advances in Bethe Ansatz 
techniques show promise for full calculations of the psi- 
nons' contribution to S(q,oj) 0,01- A direct compari- 
son of such calculations to our data would clearly be very 
interesting. 
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